o 
o 

(N 



Quark- antiquark composite systems: 
the Bethe-Salpeter equation in the 
spectral-integration technique in case of 
different quark masses 



A.V. Anisovich, V.V. Anisovich, V.N. Markov, 
■ M.A. Matveev and A.V. Sarantsev 



oo : February 1, 2008 

. Abstract 

, The Bethe-Salpeter equations for the quark-antiquark composite 
systems with different quark masses, such as qs (with q = u,d), qQ 

■ and sQ (with Q = c,b), are written in terms of spectral integrals. 

. For mesons characterized by the mass M, spin J and radial quantum 

I number n, the equations are written for the (n, M^)-trajectories with 

Q_i' fixed J. The mixing between states with different quark spin S and 

^ • angular momentum L are also discussed. 

1 Introduction 

X 

H 

5^ i The relativistic description of quark-antiquark states is a necessary step 
for meson systematics and the search for exotic states. The standard way 
to take account of relativistic effects is to use the Bethe-Salpeter equation 
[1]. Different versions of the Bethe-Salpeter equations applied to the de- 
scription of quark-antiquark systems may be found in [4, 5, 6, 7, 8]. In 
the present paper we develop the approach suggested in [2] for the Bethe- 
Salpeter equation written for quark-antiquark systems in terms of spectral 
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integrals. In [2], the systems of quarks with equal masses have been consid- 
ered such as nn, ud^ dd, ss. In this paper, the systems with unequal masses 
like qs {q = u, d) and qQ with Q = c,b are treated. 

Detailed presentation of the spectral integration method as well as the 
emphasis on its advantages was given in [2], so we need not to repeat our- 
selves. Let us only stress a particular feature of the method: it is rather 

easy to control the quark-gluonium content of the composite system that 
is rather important for the search for exotics. Another advantage consists 
in the easy treatment of the systems with high spins. 

Similarly to [2], we present here the equations for the group of states lay- 
ing on the (n, M^) trajectory. Such trajectories are linear, they are suitable 
for the reconstruction of interaction between quarks at large distances. We 
hope, by investigating high-spin quark-antiquark states, to obtain decisive 
information on the structure of forces in the region of r ~ Rconfinement- 

In this present paper we present final formulae for the {n, M^)-trajectories, 
the details of the calculations may be found in [2]. It should be immedi- 
ately emphasized that the case of different masses requires more cumber- 
some calculations. In particular, the mixing of states with {J — L; S — 0) 
and {J — L;S — 1) should be accounted for (here L and S are the orbital 
moment and spin of quarks, respectively). 

Note that for systems with the equal quark masses we have already 
obtained numerical results for a set of the (n, M^)-trajcctories such as: 
ao, ai, 02, 71", p, bi. So we hope that we have elaborated rather efficient tech- 
nique allowing us to find out realistic wave functions for the quark-antiquark 
systems. Certain aspects of numerical solutions of the Bethe-Salpeter equa- 
tions in terms of the spectral integrals are discussed in [9] . 

The paper is organized as follows. In Section 2 we define the quantities 
entering the Bethe-Salpeter equation: for equal masses they were intro- 
duced in [2], now we expand the definition for unequal masses. In Section 
3, the equations for the (n, M^)-trajcctorics arc considered for two different 
cases, for ( J = L; 5* = 0) and {J = L; S = 1) states. The technicalities 
related to the trace calculations of loop diagrams as well as trace factor 
convolutions are considered in Appendices A and B. 
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2 Quark-antiquark composite systems 



In the spectral integral technique, the Bethe-Salpeter equation for the wave 
function of the qiq2 system with the total momentum J, angular momentum 

— * — * 

L — \ J — S\ and quark-antiquark spin S can be schematically written as: 



/ 



where the quarks are mass-on-shell: kf — k'l — ml and k^ = k'2 = ml. We 
use the following notations: 

k^^{ki-k2), P^ki + k2, k' ^^{k[-k'2), P'^k[ + k'2, (2) 



P — s, P — s , — fcj/gfj^^, k^ — kj^g 
The phase space integral is written as 



ds'^^d^^\P' -k[- k'^) = f%$(s'), (3) 



2s' 



gl2 _ _ ^2^2 

The wave function reads: 



Here Qlf^'.bif]{k±) are the moment-operators for fermion-antifermion systems 
[3] defined as foUws: 

«i.(^x)=^754^1,,(A:^), (5) 

Qi;'r,Y\k^)-i^zlf-ljk^), 
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where 



k 



(2J-1)!! 
J! 



h-'- h-'- h-'- h-'- h-'- — 

/Ul H2 H3 H4 MJ 



2J 



^1 



(2J- l)(2J-3) 

+fl'/il/i2 5'/i3;U8^;U4^M6 ■ ■ ■ ^/ij + ■■■)+■■■ 



2J- 1 



L2 



\i=l 
J 



1) 



(6) 



The potential operator can be represented as a sum of the t-channel oper 
ators: 

V {s, s', {k^k'j) = J2 v}"^ {s, s', {k^k'^)) di®di, 



(7) 



Oi = I, 7^, i(7^„, ilixlb, lb- 

To write the spectral integral equations we are to transform the t-channel 
potential operator V (s, s' , {k±k'j_)) into the s-channel ones as follows: 

V {s, s', (k^k'J) = E E ik±k'±)) Ci, (Fe Fe), (8) 



/ c 



where C/c are coefficients of the Fierz matrix: 



/II 1 

'44 8 



c 



1 

3 



1 1 1 



1 i ^ 
I -1 



-i -1 



1 1 



4 8 4 4 



(9) 
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Here the summation is assumed in the ia^i^ ia^i^ structure for all indices. 



Denoting 



K {s, s', {k^k'^)) = (^^-^D) , 



(10) 



we have 



V (s, s', (A;^A;^)) = ^(F, ® K (s, (A;±A;^)) 



(11) 



= (I ® I) (s, s', (A:±A;x)) + (7^ ® 7m) (s, s', (/cx/c^)) + (ia^ 

x1/t (s, s', (/c±A;^)) + (i7M7/.®^7M7!.) (s, s', (/c±/i;i)) + (75^75) Vp (s, s', (/c±/i:i)) 

Let us multiply Eq. (1) by the operator Q^'^'.^^'j'(fc±) and convolute over the 
spin-momentum indices: 



(S,L,J) 
(n)/ii.../ij' 



(s - M^) Sp 



Y: Sp [Fe (^1 + mi)gi^'.^i)(fcx)(-^2 + m2) 



+ mOgiti'i (^■±)(-^2 + m2)l (12) 



(27r) 



$(s') X 



X V; is, s', (k^k'J) Sp [{k[ + mi)Fe (-fc^ + m2)§[y^';^)^^(^;) 

We have four states with the qiq2 spins S — and S — 1: 

1) 5 = 0; L = J, 

2) ,5 = 1; L = J + 1, J, J - 1, 

which are mixed and form two final states. The wave functions read: 



for S^0,1,J^L, 

where Cj and are the mixing coefficients with i — 1,2. 
for5 = l, L = J±1,J, 

(1,J-1,J) /I N , R 



(13) 



(n)/ii...Atj 



(«)M1---MJ 



(14) 



where and 5j are the mixing coefficients with j = 1,2. 
These wave functions are normalized: 



d^k± 
(2^ 



^{S',L' J') 
(n')Mi---/iJ 



{kj_)(ki + mi)^ 



!f&'L(^^)(-^2 + m,) 



(15) 



-l)'^^5'5^L',Lj5j'j5n'n ■ 
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3 Equations for (n, M^) trajectories 

In this section we write the trajectories for [J = L, S = 0, 1) and (J = 
L ± 1, 5* = 1) states. 



3.1 The equation for the (5* = 0, 1, J = L)-state 



There are two equations for the two states with 5" = 0, 1 and J — L. Their 
wave functions are denoted as Cj^lV'/^ „ (k±) + Di^9^'^'f^ „ (/cj,), with i = 
1,2. These wave functions are orthogonal to each other. Normahzation and 
orthoganahty conditions give three constraints for four mixing parameters 
Cj and Dj. 



Each wave function obeys two equations: 



{s - M') Xlfl^^ik^) Sp 



X 



c 



(16) 



X 



X 



(s, s', {k^k'^)) Sp 7a'(A:'i + mi)F, (-A;^ + ma) 



X 



and 

(s - M^) Sp,,,,,,P,,Zi-^J^^...^_^^,^{kj_) Sp [i75 7a (^1 + mi)^^{-k2 + ms)] x(17) 
^Si-Mj,<^3(^^) E -^^^ [^c (Ail + mi)7^,(-^2 + m2 



X y ^s')Vc{s,s',{k^k'^))Sp[t-f,^i,{k[ + m,)F,{-k'^ + m2) 



{2n) 



X 
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Now consider the left-hand side of the equation. Using the traces written 
in Appendix A and convolution of operators from Appendix B, we have: 



Ki^...^.j{k±)Sp h^{ki + mi)i75(-fc2 + 7712) 



^avii>2UzPuiZ[,2i^ ^^^^^{ki_) — 0(18) 



Here and below we use the following notations: 6 = m2 — mi, a — m2 + mi. 
Also the left-hand side of (17) contains two convolutions: 

e(3u,u,usP.^Zl^X^^^^,u,{k±)Sp [i75(^i +mi)7^(-^2 + m2)] Xj/^l^^{kj_) = 0,(19) 

^|3^^lU2U3Pu,Zi■^l,...^J,^.s{k±)Sp [jaih + mi)jj{-k2 + TTls) 

><ea...2uzP.AiL.u...Sk±) = -2s{s - A^) y,ty a(J)fcf 



X 



The right-hand side of the equation is calculated in two steps. First, we 
summarize over c: 



A{s,s',{k^k'^))= Yl A,{s,s',{k^k'^))V,{s,s',{k^k'^)) 



(20) 



^ ^ Sp Fc{ki + mi)i'y5{-k2 + m2) Sp i'y5{k[ + mi)Fc{-k'2 + 7712) x 

c=T,A,P 

xV,{s,s',ik^k'^_)) , 

Bpw {s, s', (k^k'J) = Wp'a' {s, s', {k^k'^_)) Vc (s, s', {k^k'^)) = (21) 

c=T,A,V,S 



Sp [Fc (ki + mi)7^,(-^2 + m2)] Sp [r^,(k[ + mi)F^ {-k'^ + 1112) 

c=T,A,V,S 



X 



xV,{s,s\{k^k'^)), 



and 



C^' (s, s\ {k^k'^)) = Yl Ca' (5, s\ ik^k'J) K is, s', {k^k'^)) = (22) 

c=T,A 
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Sp Fc{ki + mi)ij5{-k2 + 1712) Sp 7a/(A;i + mi)Fc (-/ca + ?7i2) 



c=T,A 



X 



In Appendix A the trace calculations are presented, and the values Ac {s, s', {k±k'j_)), 
Cc (s, s\ {k±k'j_)) are given. So after the summation, Ac, Be, Cc are written 
as follows : 

A{s,s',{k^k'^))= Yl Ac{s,s',{k^k'^))Vc{s,s',{k^k'^))= (23) 

c=T,A,P 



= -16{k^k'j) \2V7^Vt {s, s', {k±k'j)) + Va {s, s', {k±k'j)) 



-4(s- A2)(s'- A2^ 



Va {s, s', (k^k'J) + Vp {s, s', (k^k'J) 



Bp'a' {S, S', {k^k'^)) = Y Wp'a' {S, S', {k^k'^)) K {s, s' , {k^k'^)) = (24) 

c=T,A,V,S 



- A'){s' - A') (vv {s, s', {k^k'^)) + 2-^ Vt {s, s' , {k^k'^))^ + 



+4y^y^^ {s, s', {k±k'J) + 2A2 Vt {s, s', {k±k'^ 



+ 



a' Vs {s, s', {k^k'^)) + -= Vy {s, s', {k^k'M + 



-l&k'^k^, 



+ 



+^y'^^ly'k'lzVv{s,s',{k^k'^)) 

/^Va {s, s', {k±k'^)) + 2A2 Vt {s, s', {k±k'^)) 

+16kPk^, [s' - A') Vv {s, s', {k^k'^))+lQk'^k'^ (s - A^) Vv {s, s' , {kjz'^)) , 
and 

Co^' {s, s', {k^k'^)) - Y Ca' ik±k'j_)) V {s, s', {k^k'J) = (25) 



=T,A 



= 8 [2Aea'kk'P' Va {s, s', {k±k'j)) + aSa'Pk'P' Va {s, s', {k±k'j)) + 
+AAea'Pkk' Vt {s, s', {k±k'jj) + 2aea'PkP' Vt {s, s', {k±k'j))] 
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Here we used a short notation: Ea'kk'P' = kjsk'^PI^ea'f^^u-^^cond , the con- 
volution of operators is performed by using equations of Appendix B and 
recurrent formulae for the Legendre polynomials: 



zPjiz) 



2J+1' 



2J + 1 



that allows us to represent the Bethe-Salpeter equation in terms of the 
Legendre polynomials . As a result we get: 

Xlil,Ak±) A («, s\ {k^k'^)) = -4a(J) (V^^^)' x (26) 



X 



^27TT V^V^ (2^/i^^T (s, s\ {k^k'^)) + A' Va (s, s', {k^k'^))) Pj-,i{z)+ 



+4^y^^V^ (2V^Vt {s, s', {k^k'^)) + Va {s, s', {k^k'^))) Pj^,{z) 



^ Va (s, s\ (k^k'j)) + Vp (s, s\ {k^k'^)) Pj{z)+ 
\/ss' J 



2J + 
and 



16^a(J) (y^y^ 



J+1 



(27) 



X 



X 



{jTi^'^'^ ~ >^^^+i(^)) (^'^ is, s', {k^k'^)) + 2A^Vt (s, s', (^±A;1))) . 
For the right-hand side of (17) we have: 

£^v,...3^.i^Si...m.,.3(^^)C/3'(s,/,(/c^A;;)) Xlfl^^k'^) = 



(28) 



16 



2J + 3 



^a(J) (v^^fc?) 



J+i 



X 



X 



-P;(z) - JPj+i(^)) (sA (s, s', (k^k'J) + 2AVI^Vt (s, s', (k^k'^))) 



J + 
and 

(^1) =(29) 
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^2JTi (V^14 (s, s', (fcx^D) + Vt {s, s', (fcx^l))) Pj+i(^) 

+(s - A2)(s' - A^) (vy (s, s', + 2-^ V^T (s, s', Pj{z)-r 



+^27Ti (^^^ ^^^^^^^ + ^"^^ ^^^^^))) ^•^-1^'') 



2J + 

Expanding the interaction block in the Legendre polynomials series, 

K («, s', {k^k'^)) = ^ {s, s') Pj{z) = (30) 



and integrating over angle variables in the right-hand side by taking account 
of the standard normalization condition j\dz/2 Pj{z) = 1/(2 J + 1), we 
have finally: 



X 



(s-M^) [(s-A>(°'-^"^)(s)Q 
= / ^p(.')2(-A;^)^^r'^)(.0C,x 

(mi+m2)2 
J+1 

2J + 



(31) 



-H{J+l) ifj^klk'l {2V7^V^'^'^ is,s') + A'vy^'^ is,s')) + 



+e( J) (. - A^)(.' - A^) ^ is, s') + is, s') 



-HiJ - 1) {2v^V^'-'^ (s, s') + A' Vy-'^ {s, s[ 



+ 



X 



X 



(mi+m2)2 

2 



J + 



- ^(2J - 4a - 1) e( J - 2a - 1) i-^kl^kfy^^^^'^ 



X 
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where 



X (s'A l^f -2-1) (s, s', )) + 2AV^vi'-'''-'^ {s, s', ))) - 
- Je( J + 1) {s'A t/y+') (s, s') + 2A'V^V^'+'^ (s, s'))) , 



P'j{z) = J2{2J-4a-l)Pj^2a-i{z). 



(32) 



The second equation, (17), reads: 

is - \s{s - A-)^-^^^i^^y''\s)D, 

7 ^P(^08^fci(-fc!)'^+^^r'^H.')^.x 

(mi+m2)^ 

X E(2^ - 4a - 1) i{J -2a -I) (-^^y^) -2(^+1) x 

X (sA V>y-2-i) (s, s', )) + 2Ax/^l/^''-'"-^) (s, s'))) - 
- Je( J + 1) (sA (s, s') + 2Av^y/+i (s, s'))) + 

+ / ^p(.o 2 (-fe!)^^i^'^-^^(^O^i 



(33) 



(mi+m2)^ 



fX 



X 



4e(J+ 1) ^7^^i^5 (v^vl^+^^ {s,s')+2A'V^'^'^ (.,.')) + 



+e( J) (« - A'){s' - A2) (^2-^ V^'^ {s, s') + {s, s')j - 

-4 e( J - 1) ^ (V^y/-^) (., .') + 2A2 V^-'^ (s, s')) 
The normahzation and orthogonahty conditions look as follows: 



ds 



(m2+mi)2 



TT 



Pis) 



{^^y''\kl)y 2a{J){-kir (s - A2) + (34) 
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2^ J{2J + 3f 



1, i = l,2, 



and 



ds 



(m2+mi)^ 



TT 



(v.(°'^-^)(A;l))'2a(J)(-A;i)-^ - A^) + (35) 



J{2J + 3Y 



= 0. 



3.2 The equations for the (5 = 1, J = L ± l)-states 

We have two equations for the two states with S — 1 and J = L ± 1. 
Their wave functions are denoted as Aj¥^^'-'''J^{k^) + Bj¥^^;^^':'J^{k^), 
with j = 1,2. These wave functions are orthogonal. Normahzation and 
orthogonahty conditions give three constraints for four mixing parameters 
Aj and Bj. 



Each wave function obeys two equations: 



X 



(36) 



X 



X 



d^k 



i2ny 



.1 

(s, s', {k^k'^)) Sp \j^,(k[ + m,)F,{-k'^ + m2) 



X 



and 



X 



{s - M') Sp [7a (^1 + mi)7^(-fc2 + m2) 



(37) 



(J-1) 

IJ.i...Hj,l3' 



•.'{k±)^Sp Fc(ki + mi)'yp,{-k2 + 7112) 



X 
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X 



Vc (s, s', (k^k'^)) Sp -f^,{k[ + mi)Fc{-ki, + ma) 



X 



Firstly, consider (36). In the left-hand side of (36) one has two convolutions: 

Ki':lpik^)SP Uih + m^)^^{-k, + m,)] X^^'J^M = (38) 



2a{J)k 



2(J+1) 



2J + 1 



L J + l 



s - A^) + Ak] 



Ki^Zpik±)Sp [l^Cki + m,hji-k2 + m^)] Z\^-}l^M = 8a( J)A;f 
Also the left-hand side of (37) contains two convolutions: 



^/ii.../ij,/3^ 



{k^)Sp [7a (^1 + mi)7^(-fc2 + m2)l 4f+?,Jfc±) = 8a( J)^?'^+^\(39) 



2a{J)k 



2(J-1) 



2J + 1 



J 



s - A^) + 4k] 



The right-hand side of Eqs. (36) and (37) is determined by the convo- 
lutions of the trace factor Bp/a' (s, s', {k±k'^)), see Eqs. (24) , with angular- 
momentum wave functions; corresponding formulae may be found in Ap- 
pendix B. Taking them into account one has for the right-hand side of (36): 



Ki^lAk±) Bp,a' {s, s', {k^k'J) X^:^^Ak\) = 4«( J) MM] x(40) 



X 



^Hi...Hj)3' 

"2J+ 1 



J+i 



'j±l^s - A'){s' - A') (vv{s, s', {k^k'^)) + 2-^ Vt{s, s', {k^k'^))^ + 

+4{s' - A^)kl Vv{s, s', {k^k'^)) + 4(s - A^)k'l Vv{s, s', {k^k'J)+ 



+16 



-^""^■klk'lVvis.s'.ik^k'^)) 



2J+1 



+4 



Vs{s, s', {k^k'^)) + ^= Vv{s, s', {k^k'^))+ 
vss' 



J 



J + 



- {VI?Va{s, s', {k^k'^)) + 2AVt(s, s', {k^k'^)))^ MMlPj{z)+ 
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J 



and 



j-i 



X 



' ^ + 2J+1 ^ 



k'lVy{s,s',{k^k'J)Pj+,{z) + 



1^5(s, s', {k±k'^)) + ^= Vv{s, s', {k^k'^)) 
yss' 



+ 

ss' Va{s, s', {k^k'^)) - 2A^ Vris, s', (k^k'J)] ^kl^k'lPj{z)+ 
J 



+ 



s'-A^ + 4- 



-k 



12 



2J+ 1 

And for the right side of (37): 



klVv{s,s\{k^k'^))Pj.,{z)y 



Z^'LA^^) Bp'a' (s, s', {k^k'^)) xl'^tZa'ik'J = 16a{J)k'l ( ^kl^kf ) x(42) 



j-i 



X 



s'-A' + 4^^k'' 



2J+ 1 



klVv(s,s',(k^k'^))Pj+i(z)+ 



a' Vs{s, s', {k^k'^)) + -= Vy{s, s\ {k^k'^)) 

yss' 



+ 



ss' Va{s, s', {k^k'J) - 2A^ Vt{s, s', {k^k'^))\ ^Jkl^k'lPj{z)+ 
J 



+ 



s-A^ + A 



2J+1 



ki 



k'lVv{s,s',{k^k'^))Pj.i{z) 



and 



zlf'iA''^) Bp'c is, s', (k^k'J) zlf-'l^Ak'J = 4a( J) MM) x(43) 



j-i 



X 



IQ^j^klk'l Vvis, s', (k^k'J) Pj+,iz)+ 



+4 



a'Vsis, s', {k^k'^)) + ^= Vvis, s', ik^k'^)) + 
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+ 



f {V^Va{s, s', {k^k'j_)) + 2AVt(s, s', (A;xA;1)))] ^^Pj{z)+ 

J {s- A'){s' - A') (vv{s, s', {k^k'^)) + 2-^ Vt{s, s', {k^k'J)) + 



J 

2J+1 



+4{s' - A')ki Vv{s, s', {k^k'^)) + 4(s - A')k'l_ Vv{s, s', {k^k'J)+ 



J 



+^^i^j—^klk'l Vv{s, s', ik^k'^)) 



In the right-hand sides of eqs. (36) and (37), we expand the interaction 
block in the Legendre polinomials series and integrate over angle variables 
j\dz/2. As a result, Eq. (36) reads: 



{s - M") 



#^^-^'^)(5)A,+ 

ds' 



V J+ 1 ' 



s - A^) + ikl) jjj^y+^^^\s)Bj 



=(44) 



/ ^pis')8{-k'iy-'i;(y-''-^\s')A, 

J TT 



X 



(mi+m2)^ 



X 



e( J + 1) k'l (. - + 4^A:i) Vy^'\s, s')- 

a^A^ 



/2 



iJ), 



-^/ss'V^A\s. s') - 2A^ ^^'^(s, s') 



ss' 



+ 



+ 



s' - A^ + 4 



/ ^p(.')2(-A;2)^+V^-^+^'^)(A;2)i?,x 

J TT 



J 



-k 



/2 



2J+ 1 



(mi+m2)2 

X [^(.-A^)(.'-A^) (v^r^) (-,-') + 2^ V^^^^^ (.,.')) + 

+4(/ - A^)kl Vy^'\s, s') + Ais - A')k'l Vy^'\s, s') + 
+ 16- 



-^fik'lvr'\s,s') 
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-4e(J) 
J 



y ss' 



+ 



J 



+ 



2J + 1 

The second equation, (37), reads: 



/2J+1 



J TT 

(mi+m2)2 



-4^(J)^lA;? 

J+ 1 



+ 



J 



- A^)(.' - A^) (v^'-'\s, s') + 2-^ W-'\s, s') 



ss' 



+16 



2J+1 



+ 



J TT 



(mi+m2)2 

xe(J+l)A;l 
-e( J) kl 



s'-A^ + A^^k" 



2J+1 
V ss' 
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-V7^vi^\s, s') - 2A^ V^'^\s, s') 



+ 



s - + 4- 



J 



-ki 



V 



2J+1 

Normalization and orthogonality conditions are : 

ds 



(m2+mi)2 

x(^(«-A^) + 4^i) + 



(46) 



and 

ds 



J+1 



s-A')+ 4ki 



1 J 



(m2— 



TT 



-P{s) 



A,A2{^^J^y-'''\kl)) 2a{J){-klY'-'^x (47) 
2J + 1 



X 



J 



fs- A^) + 4A;1 + 



+(^^2 + ^Si)V'^^-''^n^l)V'i''^+''^H^i)8«( j)(-^i)^-'+'^+ 

+Sii32 (V^i^'^+^"')(^i))'2«(J)(-A;i)(''+^) - A^) +4;^^ 



J + 1 



0. 



Let us emphasize again: all the above equations are written for J > 0. 



4 Conclusion 

We have presented the Bethe-Salpeter equations for the quark-antiquark 
systems when the quark and antiquark have different masses. The main 
difference from the equal mass case is that there is the mixture of states 
{J — L, S — 0) and {J — L, S — 1) and that is proportional to the 
quark mass difference. The mixing between {S — 1, J — L) and (5" = 0, 
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J = L) states give rise to a strongly correlated system of equations. In the 
equation for states with total spin J we need to know all lower projections 
of the potential on the Legendre polynomials, not only the J + 1, J , J — 1 
ones. The numerical study of this equations is now in progress. 

We are grateful to L.G. Dakhno, M.N. Kobrinsky, Y.S. Kalashnikova, 
B.Ch. Metsch, V.A. Nikonov , H.R. Petry and V.V Vereshagin for stimu- 
lating and useful discussions. The paper is supported by the RFBR grant 
no 01-02-17861. One of us (V.N.M.) is supported in part by INTAS call 
2000 project 587 , RFBR 1-02-17152 and "Dynasty Foundation". 

5 Appendix A: Traces for loop diagrams 

Here we present the traces used in the calculation of loop diagrams. Recall 
that in the spectral integral representation, there is no energy conservation, 
s 7^ s', where = s, P'^ = s', but all constituents are mass-on-shell: 

l2 _ 2 i,2 _ 2 7,/2 _ 2 1,12 _ 2 

We have used notations for the quark momenta: 

K^\{k,-k2U K-\{K-k'2)., (48) 

/,-L _ u ± lJ± _ U' -L 

_ ml-ml ^ _ ml-ml 

and for the quarks masses : 

A = ni2 — mi, a = m2 + mi, . (49) 
We work with the following definition of the matrices: 

5.1 Traces for the 5 = states 

For the S — Q states we have the following non-zero traces: 

T'p = Tr [i75(^i + mi)75(-^2 + ^2)] = 2i(s' - (ma - mif), 
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i-1^{k[ + mi)i7^75(-A;2 + ma) 
T!^ = Tr i^5{k[ + mi)i(T + ma) 



2k'^{m2 - mi) + P^(m2 + mi) 



and 

Tp^Tr 
TT = Tr 



«75(-^2 + "^2)75(^1 + "^i) = 2i(s - (m2 - mi)^), (51) 

hsi-h + m2)i7^75(A;i + mi) = 2 [2A;^(m2 - mi) + P^(m2 + mi)] , 
ijbi-h + rn2)ia^^{ki + mi) 



Aie 



Ixi/apPak/S- 



The convolutions of the traces Ap — [TpT'p), A a 
(TtT^) are equal to: 



Ap= -4{s - A^){s' - 

Aa = -leA^y^z - 4-^(s - A2)(s' - A^), 

At = 



{TaT'^), At = 
(52) 



-?,2\/ss'\lk\k'lz. 



5.2 Traces for the 5 = 1 states 



For the [S = l)-states, the traces are equal to: 

T's = Tr [ii,{k[ + mi)(-A;^ + m2)] = 4A;;,i(m2 + mi), (53) 
Tv = Tr [7i(^i + mi)7^(-^^ + m2)] = 2 (g^,^{s' - {m^ - mi)^) + Ak'^k'^ , 
T'a = Tr Yi^,{k[ + mi)i7^75(-^^ + m2)] = ^.apk'^P'p, 
'li,{k'i + mi)i(T^i,{-k2 + 1712) = 
2(m2 - mi) (^i^/c;, - g^,J^'^) + (mi + ma) (g^^^Pl - ga'^^P'u)] > 



= Tr 



2i 



and 



T5 = Tr 
Ty^Tr 
TA^Tr 



l^'{-k2 + m2){ki + mi) = 4fc^,(mi + m2), 



(54) 



lp'{-k2 + m2)7^(A;i + mi)] = 2 [5-^^,(5 - (m2 - mi)^) + 4A;^,A;^ 
7^,(-A;2 + m2)niu^5{ki + mi) = -^epf^apkaPn, 
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Tr^Tr 



Ip'i-h + m2)ia^^{ki + mi) 



= 2i 



2(m2 - mi){gjr^K - gj'M + ("^2 + mi) (gj,^P^ - g^f^P^) 



Corresponding convolution Bs = {TcT^) reads : 



{Bv)p'a' - 4 {gj,^,{s - A^)(/ - A^) + 16 kj.k'^Jklkf z+ 



.2 A 2 



+4 k^,kt,,{s' - A^) + 4 k'^k'^{s - A^) + 4 k^,k\ 



{BA)p'a' = -IQVSS' 



{BT)/3'a' — —8 



(55) 



6 Appendix B: Convolutions of trace factors 

Here we present the convolutions of the angular-momentum factors. Let z 
be: 



jk^k'^) 

_L V _L 



Analogous convolutions for (5" = l)-states are written as follows: 



a{J) 

J+1 



klJk'?] X 



X 



'^Ap,,,^,(^) kjk^ + ^^Bp,^^^,{z) k'p^k'^+ 



+Cp,,,^,(^) + Dp^^j^z) k'p^kt + ( ^Jkl^k'l ) Ep^^^^,{z) g^^ 



'-i-i 



(56) 



The convolutions for the {S = 1) states read: 



(58) 
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X 



X 



(59) 



fV I A / fV I 



J-2 



X 



(Ap,,,^,(^)-(2J+1)^^(^)) kjki+ 



{Bp,^,^,{z)-{2J + 1)Aj{z)) 



+ 



+ fcp.,.,,(^) + ^^^^^Pj(^) - 2(2 J + l)Bjiz)] k^k'^ + 



J+1 



+L>P,,,^,(z) + {^kl^k'l)Ep^^,^,{z) gj, 



and 



-a(J) (Jj^Jk^y \PP') X 



(2J + 3)2 / ^ ^xJ-i 



(^ + 1)= 



X 



k'? , , a/A-;^ _^ _^ ^ ^ 



A/ I \ / A/ I 



" -M J ( ) ^"'^l ''2^3 ^l^ijll • ■ -/i J ,1-3 ( ^± ) 



(61) 



(62) 
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J -I 



Here 



2^Pj(;^) + [J^2_(j_^2)]Pj+i(;^) 



(63) 



(1-^2)2 

[(1 - -/)^^ + (-/ + 1)] PA^) + [(2-/ + 1).^ - (2.7 + 3)] ,P,^,{z) 

(1-^2)2 

[{J + 2)z''-J]Pj{z)-2zPj+i{z) 

(1 - ^2)2 
^Pj(^) -Pj+l(^) 



(1 - ^2) 



For the factors 
.(j+i) 



■MJ/3(^^)^S2^-- MJ"(^^)^"' ^^^^^ = A;, k', we need more compli- 
cated convolutions, namely: : 



J+1 



a{J) JklJk'l 



J+2 



2J + 1 
J+1 



zPj+i{z) 



J 



J + 1 



2J+ 1 



J+1 



as well as for the factors 
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^^^iS™/3(^^)^i(;2...,a(^l)^a = kla{J) y^i^l^' Pj{z) , 

and for the factors 



(j-i) 



/3 ^,,,..:.,,,/3(^±)4f;5.;.,,a(fcl) k'^ = a{J) [^Jklsjk'l) Pj{z) , 

ik±)Zl'-Xj,aik'±) k'^ = k'la{J) (^^^y~' Pj-,{z) , 



k'(3 zliiJ.2-tJ.jAk-i-)ZiitJ.2-tJ.jAk'±) ka 



J2 



2J + 1 



J 



= a{J) [\JklsJk 



zPj-i{z) 



J+1 
J 



and for the factors Kpep^^^^^^P^,Zi/^l^...^j^^^{k±) eaAiA2A3^{i4i...Mj,A3(^J 



^/3i^li^2i'3-^l'l-^f'2/tl---MJ,I'3( 



I/l 1/2 1-3 -fl/l ■^i'2Ml ■ ■ -M J -i^a 

(2J + 3)2 



-(J) 

A2Ati- 


•^j,A3(^±) = 


0, 


-(J) 
A2Ati- 


•Mj,A3(^±) k'a = 


0, 


-(J) 
A2Ati- 


■IJ.j,\3^k'±) k'^ = 


0, 


A2/il- 


■IJ.j,\3^k'±) ka 





-a{J) 



rv I \/ 



,7+1 



(PPO [;.Pj(;.)-Pj+i(;.)] , 



(J+l)3 

9pa^Pi'ii'2i'3Pfi'Zii]ii-lxj,V3{k±) ^aXi\2\3P\iZ\2Hi-HjM^k'l) — 

J(2J + 3)2 



a(J) JklJk'l {PP')Pj{z) . 
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